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" Abstract 

We consider diffusion of a cold Fermi gas in the presence of a random optical speckle potential. 
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The evolution of the initial atomic cloud in space and time is discussed. Analytical and numerical 
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results are presented in various regimes. Diffusion of a Bose- Einstein condensate is also briefly 
discussed and similarity with the Fermi gas case is pointed out. 
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I. INTRODUCTION 



Transport of cold atomic gases in the presence of a quenched random potential is a rapidly 
developing field of research [lj. In a typical set-up the gas is released from a harmonic trap 
and undergoes expansion, while being scattered by the random potential. At some later 
time an image of the expanded atomic cloud is taken and, thus, information about the mode 
of transport (ballistic, diffusive or localized) is obtained. The random potential for atoms 
is obtained by creating a random pattern of light intensity (optical speckle). Experiments 
on propagation of cold atoms through optical speckles have been limited so far to one- 
dimensional (Id) geometry and have culminated in observation of Id Anderson localization 
for a Bose-Einstein condensate (BEC) [2J[3]. 

There is a considerable amount of theoretical work on diffusion and possible localization 
of an expanding BEC cloud in two and three dimensions j3HS] - The same problem can be 
also addressed for a cold Fermi gas - a system which is intensively studied in recent years (see 
[TU| ITT] for recent reviews). Diffusion of an expanding Fermi gas, in the long time limit and 
for a Gaussian white noise potential, was discussed in [12j. In the present paper we consider 
the experimentally relevant case of a speckle potential, concentrating on 2d geometry. In 
Sec. II we write down the basic equations which govern the evolution of a diffusing Fermi 
cloud. In Sec. 777 we summarize, following [13], the behavior of the diffusion coefficient 
D(k), as a function of the particle wave number k, in a 2c? speckle potential. In Sec. IV 
we study the density n(f,t) of a diffusing Fermi g function of position and time. 

Since n(f,t) is expressed by an integral which cannot be calculated analytically, we resort 
to numerics in combination with an analytic treatment of some limiting cases. In Sec. V we 
briefly discuss the evolution of the shape of a diffusing BEC and point out some similarities 
(and differences) with the case of the Fermi gas. 

II. BASIC EQUATIONS 

We consider N fermions at zero temperature, initially trapped in a harmonic potential. 
At time t = the trap is switched off, while a random potential V(r) is switched on. Our 
aim is to study the dynamics of the atoms, upon their release from the trap, in the presence 
of the random potential. In many circumstances interactions between the fermions have only 
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a minor effect on their dynamics. This is particularly true for a polarized Fermi gas when 
the Pauli principle eliminates the main mechanism (the s-scattering) for the interaction. In 
the absence of interactions the single particle wave functions, \I/ n (r, £), describing individual 
atoms, evolve according to: 

ihd t V n (r,t) = -^-A* n (r,t) + V(r)S> n (r,t) , (1) 

with the initial condition \l/ n (r, 0) = $> n (f), where $ n (j) is n'th eigenstate of the harmonic 
potential ^mu 2 r 2 , and V(r) is the random potential, with zero mean and a two-point cor- 
relation function (V (f) V (r^)) = T (r*2 — fx). 
The formal solution of Q is 

* n (r,t) = j dRG(f,R,t)$ n (R) , (2) 

where G(r,R,t) is the retarded Green's function of the Schrodinger equation ([!]). The 
quantum expectation value of the particle density (per one spin component) at time t and 
for a given realization of randomness is 

< n(f,t) >=J2fn\*n(ft)\ 2 = j dRdR'G*(f, R, t)G(f, R', t) ^ fn§* n (R)$n(R') , (3) 

n n 

where f n is the occupation function, which for zero temperature is given by the step function 
Q(Ep — E n ). Averaging < n(f,t) > over the disorder yields 

< fi(f,t) >= f dRdR'G*(f,R,t)G(f,R',t)J2fnK(^n(R') . (4) 

J n 

In order to average the product of the two Green's functions in Q we first Fourier transform 
to the energy representation 

de f dQ int „ l„s„,^ ^» 1 



G*(fR,t)G(r,R>,t)= - —e~G*(r,R,e + -Q)G(r,R>,e--n). (5) 

The product in r.h.s. of (|5]), in the diffusion approximation, is represented diagrammatically 
in Fig.l. 

R 
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Figure 1: Diagrammatic representation of the product in Eq.(5) 
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The two straight lines represent a wave and its complex conjugate propagating from their 
sources (R' and R) to an intermediate point fi. At this point the two waves "recombine" and 
the wave intensity propagates by diffusion to the observation point r. The corresponding 
algebraic expression is 



G*(f,R,e + ln)G(f,R',e-ltt) = - f df 1 P £ (f,rl,fi)G*(f 1 - J R,£+^)G(f 1 - J R', £ -^), 

2 I T £ J 1 1 

(6) 



where r £ is the scattering mean free time at energy e, P e (r, r*i, fi) is the diffusion ladder [TJ 
and G(f, e) is the average Green's function. For the latter Q can be neglected, in comparison 
with e, and its explicit expression is 



G(f, e ± -n) « G(f, e) = G (f, e)e" 



(7) 



where Go is the free Green's function and l e is the single particle mean free path. Since the 
Green's functions in Eq.(|6]) rapidly decay (at a distance l £ ), the slow varying diffusion ladder 
P £ (f,fi,Q) can be taken out of the integral, with the argument f\ being replaced by R ~^ R . 
Performing the remaining integral and returning to ^ yields 



G*(f, R, t)G(f, R',t) = -- f deP £ I f, R + R ,t j ImG ( /? - /?'. c 



(8) 



where the diffusion propagator 

P £ (f,R,t) 



exp 



f- R 



\ 



AD £ t 



(9) 



/ 



{^D £ tf 2 

is the Fourier transform of P £ (f,R,Q) and D £ is the diffusion coefficient at energy e. The 



necessary condition for the above derivation is kl £ ^> 1, where fc = a/ 2me / h 2 . 
Substituting |8) into (4) and using the fact that for weak disorder —^ImG(k,£) ~ 5(e — e^) 
one obtains 



< n{f, t) >= n{f, t) = j dR j dpP p {f, R, t) ^ f n W n {p, R) 



(10) 



where P p is given by ( 9 ) with e = £— and 

1 



w n (p,R) = J dpe^K{R+\m n {R-\p) 



(11) 
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is the Wigner transform of <E> n (f) . In the classical limit (n 3> 1) the Wigner function 
W n (p, R) for an eigenstate n becomes |15j 

1 



(12) 



where E n is the energy of state n and v{E n ) is the density of states for a particle in a 



harmonic trap. Substituting this into (10) and replacing summation over n by integration 



over energy up to the Fermi energy Ep we finally obtain 



n(f, t) 



dR 



P 



2 1 



— J —PJf,R,t)Q(E F 



muj 2 R 2 ). 



(13) 



III. DIFFUSION COEFFICIENT IN SPECKLE DISORDER 



In order to proceed with the evaluation of the integral in (13), an explicit expression is 
required for the diffusion coefficient 

hkl B 



D(k) 



dm 



(14) 



where Is is the Boltzmann transport mean free path and k = f. So far we have not 
specified the type of disorder. Now we specialize to a two-dimensional (d = 2) speckle 
potential, generated by transmitting laser light through circular diffusive plate, whose two- 
point correlation function is given by j5] 

2 



r(fi - r 2 ) 



4V 2 



Ji(k 


In. - f 2 |) 


ko | 


Fi - r 2 





(15) 



where J\ is the first-order Bessel function, Vo is the standard deviation and ko is the inverse 
correlation length of the random potential. The latter is related to the laser wavelength and 
numerical aperture of the imaging device. Then, in the weak disorder limit, the mean free 
path is given by |13j 

„ 2tt 



1 

kl B 



—dOT ( 2^ 
2tt V k o 



sin 



cos(0)) 



(16) 



where 7] = |f is the measure of the potential fluctuations strength, E = is the "corre- 
lation" energy and 



r (k) 



i 



8 arccos , 

V V2 



K 



e (2 - k) 



(17) 



In the limiting cases, k <C ko and k ^> ko, (16) may be approximated as [1 3J : 

I 

kl B ~ 



Airvp \ko J 

45?r / jA 5 
128»? 2 \k J 



k < k 
k > k 



(18) 



In Fig. 2(a) we compare approximations (18) to the exact numerical evaluation of (16). 
The optimal choice of a point, separating between the two asymptotics, is the crossing point 

1 /3 

k cr = Xko, with A = (1P2) ~ 0.41. With this choice (16) is approximated as: 

kl B = 



45tt / _fc_ 
128r7 2 \k 



k < Xk 
k > Xk 



(19) 



The approximation (19) differs from the exact numerical solution of (16) by a numerical 
factor of order unity. This is demonstrated in Fig. 2(b), which shows the ratio between the 
two. 
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Figure 2: (Color online) (a) Comparison between the numerical (o), Eq, 16, and approximate 
(lines), Eqfl~9~| solution for the Boltzmann transport mean free path The small k (solid 
line) and large k (dashed line) asymptotics cross at k cr = Xk^; (b) The ratio between 
the exact numerical solution of ( Jl6| and the approximate expression ( fj~9| ). The maximal 
deviation, obtained at the crossing point k cr = Xk , is about 1.76. 
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Thus, we write the diffusion coefficient as 

[D < (k) = ^(f] 2 , k c <k<\k 
D(k) = { Zr^s , (20) 

{D>(k) = ^{±) , k>Xk 

where Dn = , ° 2 . Note that we have introduced a lower cutoff k c , which is derived from the 
Ioffe-Regel criterion k c ls = 1 (see below). For k < k c the diffusion approximation employed 
in this paper is not valid any more: thus, particles with k < k c remain localized in the 
vicinity of the initial trap. 



IV. EVOLUTION OF THE DENSITY IN SPACE AND TIME 



Using the above explicit expression for the diffusion coefficient, one can calculate the 



atomic density profile n(f,t) (13). It is convenient to introduce the following dimensionless 
variables: 

/ 

r = p J — 

tt-M ax 

t= h , (21) 



< 



n 



nR 2 



Max 



where R Max 



2Ep 



/8N 



is the initial size of the atomic cloud and to 



Max _ 2 V 2 V2N 



IS 



a characteristic diffusion time. Let us note that for t — > the diffusion kernel ^ becomes a 
delta function 5(f — R) and the density approaches its initial shape of the inverted parabola, 

2N 



n (r,0) 



7T 



-(l-? 2 ) 



(22) 



which corresponds to the Thomas-Fermi approximation for N fermions in the harmonic 



trap. Since it is difficult to calculate analytically the integral in the expression (13), below 
we consider various special cases. 

In the long time limit the atomic cloud will spread to a distance much larger than its 



initial size Ruax- Then, one can set R = in the diffusion kernel in (13) and integrate over 
R, with the following result: 



where 



n(f, t) 



$(p) 



(27rft) : 



;P p (f,0,t) $(p) 



7TR 



Max 



g ) Q (pf - p) 



(23) 



(24) 



is the momentum distribution of the gas. Eq. (23) has a simple interpretation: it describes 



classical diffusion of particles with momentum p and energy e = 2— and with a momentum 



dependent diffusion coefficient given in (20). It is interesting to note that (23) is completely 
analogous to the corresponding expression for a diffusing BEC, with kp = Pf/H being 
replaced by the inverse healing length l/£ |H E]. 



The integral in (23) cannot be calculated analytically due to the complicated dependence 



of the diffusion kernel on the particle momentum p. A considerable simplification occurs if 
one assumes pp <C hk . In this case all atomic wave numbers satisfy the condition k < k so 
that correlations in the random potential do not come into play. The diffusion coefficient is 



given by D < (k) (see (20)) in the whole range of integration which corresponds to the limit 



of an uncorrelated, white-noise potential. The expression (23) reduces to: 

" kF kdk 



n(r, t) 



R 2 



Max 



8nt 



D<(k) 



cxp 



4D<(fc)i 



kp 



(25) 



Let us stress that the white noise limit, Eq. (25), requires that the typical strength Vq of the 



random potential must be smaller than the correlation energy Eq, so that the parameter t] = 
|^ <C 1 [5j. Indeed, the white noise condition, k <C ko, is compatible with the weak disorder 



requirement, Mb > 1, only if r\ -C 1 (see (19)). This inequality implies k c = V Anrfk^ <C k { 



Furthermore, in order for the weak disorder requirement to be satisfied for the great majority 

v 2 

of the fermions, we must require kp ^> k c , i.e. Ep 3> Switching to the dimensionless 
variables and performing the integral yields: 
2Ns 



n(r, t) 



t 



1 £_ 

+ 2-Ksrfe 2 * + I 1 + 



~-2 



Ex 



irsr 



2r] 2 t 



, (26) 



where s — and the special function E\ (x) is the exponential integral |16| . The afore- 
mentioned condition Ep ^> ^ implies that the parameter srj 2 < 1. As an experimentally 
relevant example, we consider the Li 6 atoms in the isotropic trap with the harmonic confine- 
ment frequency ^ ~ 160 Hz and the speckle scale |^ = 0.5/im. For r\ = 0.05 and s = 12, 
this corresponds to iV ~ 10 4 atoms trapped in the initial cloud of the radius R max ~ 50fim 
and the typical time to ~ 0.7ms, which is about two orders of magnitude larger than the 



Boltzmann transport mean free time tb- Expression (26) is plotted in Fig. 3 for s and 



r\ specified above. In Fig. 3(a) n(r,t)/N is shown as a function of normalized time and 
distance. The chopped part of the plot corresponds to the region where the approximation 
of long time limit is not valid. Fig. 3(b) depicts snapshots of the density at different times. 
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Figure 3: (Color online) (a) Normalized atomic density n(r,t)/N Eq.(26), as a function 
of normalized time and distance for s — 12 and rj = 0.05. The chopped part of the plot 
corresponds to small times and r for which the approximation of long time limit is not valid; 
(b) snapshots of the density for s — 12 and r\ = 0.05 at four times ^ = 400, 500, 600, 3000. 



Let us discuss the obtained expression (26) in different regimes. For r < Rmclx (i- e - r 1) 



and for large times t>\, using the expansion of E\(x) for small values of x [TB] 



£i(x) = - In £ + 0(1) , 



(27) 



Eq.(26) simplifies to 



~ N 2Ns , 
n(r, t) ~ ~ In 



2TTsrf 



(28) 



In the main region, Kuax < r < wAD^kp^ (i.e. 1 < r < 



expression (26) is not 



intuitive and, for the visualization, in Fig. 4 we compare it with the solution for constant 



D = D < (k F )(see Eq.(20)) 



n(r, t) 



NS TTsf' 

t 



[1 - 2nsr] 



2\2 



(29) 



where the factor in the parentheses accounts for the lower momentum cutoff k c . As expected, 
the solution for 2d speckle has a more compact shape and the density decays faster than in 
the case of constant D, for which the density shape is Gaussian. 
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Figure 4: (Color online) Comparison between the Fermi gas density for 2d speckle, Eq. (26) 



and constant diffusion coefficient D, Eq. (29), at t = 400. The parameters are s = 12 and 



7] = 0.05, for which the momentum cutoff p c <C Pf- 



For larger time, such that r <C ^AD < (kp)t (i.e. r <C y -^), and for r > Rmo,x, Eq.(26) 



with the help of (27), reduces to 

n(r,t) p 



2iVs 
t 



In 



mm 



1 



TTST 2 2TIS7]' 1 



(30) 



which differs from the "usual" large time 1/t behavior by the logarithmic factor. The later 
originates from the diffusion constant dispersion. Finally, for r > ^yAD < (kp)t (and for 
r < 4I ^j ) where (26) is still valid), one can use the large x asymptotic expansion 



E 1 (x) = x- 1 e- x [l + 0(-)] 



to obtain 



n(r, t) 



2Nt 
r 4 sir 2 



exp 



~5 
ixsr 



(31) 



(32) 



which differs from the Gaussian decay by the algebraic factor =%. Let us note that this 
asymptotics is for zero temperature, i.e. when there is a sharp cutoff of the atomic momen- 
tum distribution at Ep. 

In order for the condition kp <C A; to be fulfilled, the number of atoms iV should be fairly 
small. When N increases, for a fixed frequency trap u, one arrives to the opposite regime 
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kp 3> ko- The integral (23) is then split into two parts (using (20)): 



n{r, t) 



f? 2 

^Max 

8nt 



Xk 



+ 



Afco 



kdk 
D^k) 



exp 



kdk 

..2 



4D>(ife)t 



exp 
1 



/c 2 



1 - 



+ 



n<(r, i) + n>(r,i) 



(33) 



The first part, n < (r,t), describes contribution of "slow" particles which diffuse with the 
coefficient D < (k) ~ A; 2 , as in a white noise potential. The second part, n>(r, t), corresponds 
to "fast" particles for which correlations in the random potential lead to a sharp increase in 



the diffusion coefficient, D > (k) ~ k 5 . For an arbitrary r, the solution of (33) is given by : 



n(r, t) 



2Ns ( (2ix7 2 \ ( 64 r 2 



t 



t 



457T t 



(34) 



where 

Fx (x) = 2irsr] 2 e~^ 
64 



X 2 s 



F 2 x) 



225vr 2 



s l , 



-e~^ + ( 1 

1 x 
5' A5 



sx 

Y 
1 

5' V2 



Ex 



x 
A 2 



Ex 



x 



5 

S\2 



X 



X 



47T?7 2 

r 



"3 


(f) § * 


-r 


3 x 




.5' 




5' A5_ 


)] 



and T (a, 2) is the incomplete Gamma function. In Fig. 5 the density n(r,t)/N is plotted 
for s = \ and r\ = 0.01 as a function of the normalized time and distance. One can observe 
that for fixed r the time evolution of the density exhibits a slight kink. It is due to the 
division of particles into two groups - "fast" (k > \ko) and "slow" {k < Xko). It is not clear 



whether this is a genuine physical effect or an artifact of the approximation (20) for D{k) 
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Figure 5: (Color online) Normalized atomic density n(r,t)/N Eq.(34) as a function of nor- 
malized time and distance for s = | and r/ = 0.01. 



For r > ^yAD > (kp)t 1 the asymptotic tail of the solution (34) can be written as 

1 



n(r, t) 



exp 



-fa' 



r 2 



t 



(35) 



where (5 = 

' 457TV2 

The above discussion pertained to the case 77 <C 1. In the opposite case, 77 ^> 1, disorder 
correlations are important for all relevant values of k, so that one should use D > (k) in the 
whole region of integration. This is because the weak disorder condition, Mb > 1, can now 



be satisfied only for k > ko, see Eq.(18). The cutoff k c , below which this condition fails, is 



now given by k c ~ k^rj^ > ko. Thus, n(r,t) is given by the second term in (33), but with 
the lower limit of integration being equal to k c . 



Let us return to the question of validity of the expression (23). It was argued that the 



transition from (13) to (26), i.e. the replacement of P p (f,R,t) by P p (f, 0,t) is justified for 



sufficiently long time. However, whether a given time t can be considered "sufficiently long" 
depends on the value of the diffusion coefficient for the relevant particles. It is clear that 



for "fast" particles, which rapidly diffuse out from the vicinity of the trap, (26) will become 



accurate at earlier times than for slow particles, which tend to stay in the vicinity of the 
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trap for much longer. Formally, the replacement of P p (f,R,t) by P p (f,0,t) requires 

fPMax 



2D(k)t 



< 1 



and 



F? 2 

11 Max 

AD(k)t 



< 1 . 



(36) 



(37) 



For r > Rmox it is sufficient to satisfy only (36), because (37) will be satisfied automatically. 
Then, for some fixed r one can identify three different time limits. For short times, t < 

r RMax 

2D(k F )f 



(36) breaks down. This, however, is of no consequence since at such small times 



even the fastest particles have not yet arrived to point r (more precisely, particle density 



there is exponentially small). For intermediate time, t 



4D(k F ) 



the fast particles arrive to 



point r and the above conditions are satisfied for these particles (these conditions are not 
satisfied for slow particles but this is irrelevant since, for these r and t, the contribution of 
slow particles to n[r,t) is small). For longer times, t > A ^ kp y the fast particles (k ~ kp) 
have already diffused away and slower particles start to arrive at point r. The arrival time 
for particles with a given value of k (smaller than hp) is of order so that the condition 



(36) is satisfied for these particles. It follows, thus, that for r ^> Rhiax the above conditions 
are satisfied for the "relevant" particles, i.e. the ones which dominate the concentration at a 
given r and t. 

For r < Rmclx the more stringent condition is (37) and in order for it to be satisfied for 
the smallest wave number k = k r , one needs t > 



R 2 

'"max 

4D(fc c ) 



i.e. t > 



(we assume here ?|<1). 



In order to obtain more accurate results for r < Rmclx and for not too long times, one has to 



return to Eq. (13) and use the kernel P p (f, R,t), rather than the long time approximation 



P p (f, 0, t). It turns out that for r = and for the case kp <C ko Eq.(13) can be evaluated 
exactly for an arbitrary time t: 



n(0,t) = 4iVs 



i] 2 I < 



1 TVS 

2t 



Ex 



TVS 



t 



_2trf 



For t > rrs, (38) can be cast in the following form: 

2Ns 



n(0,t) 



t 



In 



mm 



t 



7TS' 2HS7] 2 



0(1) 



(38) 



(39) 



It is instructive to compare the above results for n(0, t) with the solution for the constant 
diffusion coefficient D = D < (kp), which, for t > its , is approximately 

Ns 



n(0,t) 



t 



1 — 2ttsi] 



2\2 



(40) 
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In the case of the speckle disorder, the decay is slowed down by the factor 



2 In 



mm 



t l 



its ' 2-irsri 2 



reflecting slower diffusion of less energetic particles. As an illustra- 



tion, in Fig. 6 we compare these two cases for s — 12 and rj = 0.05. Note that for t > \, 



(38) reduces to Eq. (28). 
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Figure 6: (Color online) Comparison between the Eq.(39) and the corresponding expression 



for a constant diffusion coefficient Eq.(40), for s = 12 and 77 = 0.05. 



For the case kp ^> k the expression for n(0,t) is more cumbersome and involves in- 
complete Gamma functions. The main differences from the case kp <C k occurs for times 
si <^t<t^l. For such times, the density decays as 



n(0,t) ~ n, 



(41) 



For larger times, t ^> 1, the behavior of the density n(0,t) will be generally similar to the 
case kp <C ko, as discussed above. 

Finally, let us calculate the variance Ar 2 (t) = J n(f,t)r 2 df of the expanding density 



profile. Substituting n(f,t) from (13), we obtain 

Ar2 (t) = 2^2d R2max + 2dDt ' 
where D denotes average of D (p) over the momentum distribution: 

2T(d) 



(42) 



D 



Trfrf 



IpI<pf 



Vf 



{p F - P 2 y 2 D(p) 



(43) 
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Similar result was obtained in Ref.[5j for the variance of the BEC cloud expanding from the 
harmonic trap. In that case, however, the momentum distribution is given by the inverted 
parabola in both d = 2 and d = 3. 



V. DIFFUSION OF A BEC 

Previous sections were devoted to a cold Fermi gas. In this section we briefly discuss 
diffusion of a BEC expanding through an optical speckle. This problem has been addressed in 
a rather detailed and experimentally relevant paper of Miniatura et al [5J, with an emphasis 
on the limiting stationary density distribution. Here we concentrate on the earlier stages of 
the time evolution of the expanding BEC cloud. Our treatment will be within the mean field 
(Gross-Pitaevskii) approximation, when the BEC can be described by a single macroscopic 
wave function *&(r, t). The expansion occurs in two stages, when the first stage is dominated 
by the nonlinearity whereas the second stage describes a linear evolution in the presence of 
disorder |4Tl6| ITT]. Initially the condensate is prepared in a harmonic trap (frequency u) and 
its energy is dominated by interactions, i.e. by the nonlinear term in the Gross-Pitaevskii 
equation. At time t = the trap is switched off and the BEC undergoes a free (ballistic) 
expansion for a time to equal to few (-). By that time the interaction energy, stored in the 
initial wave packet, is converted into the kinetic energy of the condensate flow so that the 
interaction can be neglected. At t = to the random speckle potential is switched on and the 
BEC evolves according to the linear Schrodinger equation, with the static potential V(f). 
It has to be solved with the initial condition \P(r, to) = $(r), where $(r) is the condensate 
wave function at time t Q . Its shape is given by an inverted parabola, with superimposed 
rapid phase oscillations indicating large kinetic energy (see Eq.(23) of Ref. [5J). Measuring 
the time from the instant to, the standard treatment leads to the following expression for 
the condensate density, averaged over various realizations of V(r) (compare to (10)): 

n B {r,t) = j dR j dpP p (r,R,t)W B (p,R) , (44) 

— * 

where W B (j>, R) is the Wigner function corresponding to the wave function $(f). Let's 



compare (44) with the corresponding expression (13) for fermions. Defining the "effective 



Wigner function" of the Fermi gas as 
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we can write (13) exactly in the form as (44), with Wb replaced by Wp. The two functions 
have much in common. Integration over p and over R, respectively, shows that the spatial 
distribution and the momentum distribution for Wp(p, R) are inverted parabolas (in 2d), 
with characteristic length Ruax = 7^7 and characteristic momentum pp. But such inverted 
parabolas (with kp replaced by the inverse healing length l/£ of the BEC prior to the release 
from the trap) are well known to correspond to the condensate wave function $(f) and, thus, 
to the Wigner function Wb(p,R)- It is therefore clear that the dynamics of a BEC and of 
a Fermi cloud must be quite similar. (This similarity has been used in [3] to propose a 
single parameter scaling for BEC dynamics). For instance, in the long time limit discussed 



in Sec. IV, when R in the diffusion kernel P p (f,R,t) can be set to zero, (44) will involve 
only the momentum distribution J dRWs{p,R) and, thus, the functional form of ns{T,t) 
will be identical to that of the Fermi gas. Therefore all the results based on Eq.((23|) - such 



as those given in (25) or (33) - hold also for a BEC (with the replacement kp — > l/£). 

One should keep in mind that, in spite of having much in common, the functions Wb(p, R) 
and Wp(p,R) are not identical (indeed, two Wigner functions with the same spatial and 
momentum distributions do not necessarily coincide!). Therefore, for r < Ruax (and for not 
too long times) the shape of a BEC cloud is expected to differ significantly from that of a 



Fermi gas. Eq.(25) is not applicable in this regime and one should use the more elaborate 



Eq.(44) which is the bosonic counterpart of Eq.(13) for fermions. 



VI. CONCLUSIONS 



We have considered diffusion of a Fermi gas in the presence of a random optical speckle 
potential. The problem, although straightforward in principle, is quite involved technically 
and it differs in several respects from the standard diffusion problem encountered in con- 
densed matter physics [14j. One difference is that a broad range of particle momenta has to 
be considered, rather than a narrow interval near the Fermi momentum (as is usually the 
case for the electronic systems). Another difference is that the speckle potential has long 
range correlations. 

We have emphasized the importance of the parameter 77 = ||, where Vq and E are, 
respectively, the typical amplitude and the "correlation energy" of the potential j 18 j . For 
77 <C 1, particles with wave number k < ko do not feel correlations in the potential and diffuse 
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as in a white noise potential. For 77 ^> 1, on the other hand, correlations are important for all 
particles, regardless of their momenta. In that case an accurate estimate of the lower cutoff, 
k c , below which classical Boltzmann transport is impossible, becomes somewhat ambiguous. 

2 

Our estimate was based on the Ioffe-Regel criterion, JcJb = 1, and it leads to k c ~ kof]^ > k$. 

4 

This corresponds to a critical energy E c ~ E^rj^ which is slightly smaller than Vq. Since, 
however, E c is much above the percolation threshold E p (in two dimensions E p = 0), there 
exists a broad range of energies in which particles can propagate by classical percolation (of 
course, in 2d, and at sufficiently large distance, quantum interference will eventually take 
over and lead to localization). Such "percolating particles" were not accounted for in our 
treatment. This omission can be partially rectified by treating k c as a phenomenological 
fitting parameter whose value is determined from experiment. 

Although the paper is devoted primarily to fermions, we have discussed in the last section 
diffusion of a BEC. It turns out that, within the Gross-Pitaevskii approximation, the shape 
of a diffusing BEC cloud is remarkably similar to that of a Fermi gas. 

All kinds of localization effects have been neglected in the present paper, so that the 
weak disorder requirement, kpls 3> 1, is a necessary condition for the results to be valid. 
Finally, we have focused on the 2d case. Similar calculations can be performed also in 3d, 



starting from Eq.(13). Of course, one has to use the 3d diffusion kernel and the appropriate 
expression for the diffusion coefficient D(k) in a 3d speckle potential. 
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